Necklace— Cloverleaf Transition in Associating RNA— like Diblock Copolymers 



M.V.Tammi'2, S.K.Nechaeyi'^ 
^LPTMS, Universite Paris Sud, 91405 Orsay Cedex, France 
^Physics Department, Moscow State University 119992 Moscow, Russia 
^Landau Institute for Theoretical Physics, 117334 Moscow, Russia 

(Dated: June 20, 2006) 

We consider a Am,B„ diblock copolymer, whose links are capable of forming local reversible 
bonds with each other. We assume that the resulting structure of the bonds is RNA-like, i.e. 
topologically isomorphic to a tree. We show that, depending on the relative strengths of A-A, 
A-B and B-B contacts, such a polymer can be in one of two different states. Namely, if a self- 
cissociation is preferable (i.e., A-A and B-B bonds are comparatively stronger than A-B contacts) 
then the polymer forms a typical randomly branched cloverleaf structure. On the contrary, if 
alternating association is preferable (i.e. A-B bonds are stronger than A-A and B-B contacts) 
then the polymer tends to form a generally linear necklace structure (with, probably, some rear side 
branches and loops, which do not influence the overall characteristics of the chain). The transition 
between cloverleaf and necklace states is studied in details and it is shown that it is a 2nd order 
phase transition. 



I. INTRODUCTION 



It would not be a strong exaggeration to say that since the very beginning the statistical physics of macromolecules 
was driven mainly by biological motivations. The problems of phase rearrangements in biopolymers, such as proteins, 
DNA and RNA molecules, have stimulated an intensive development of new theoretical approaches aiming to describe 
and to predict the structural, conformational and functional changes in biopolymers under variation of external 
conditions (see, for example PjSHI)- 

A very important biological role has a wide class of so-called "associating" polymers. The associating polymers, 
besides the strong covalent interactions responsible for the very chain-like structure of a macromolecule with frozen 
sequence ("primary structure") of monomer units in the chain, are capable of forming additional weaker reversible 
temperature-dependent (i.e. "thermoreversible" ) bonds between different monomer units. Many biologically impor- 
tant macromolecules, like proteins and nucleic acids, belong to the class of associating polymers. It is known that just 
the presence of thermoreversible bonds is crucial for formation of various equilibrium spatial conformations of such 
macromolecules (so-called, ternary structures), and, hence, affects biological activity of these macromolecules Q. 

The physics of associating polymers is typical for statistical mechanics of complex systems and consists in an 
interplay between several entropic and energetic factors. In the particular case of associating polymers one can pick 
out the following three major factors, which, together with the primary structure of the chain, give rise to all the 
variety of possible thermodynamic states in these systems: 

- the energy gain due to the formation of thermoreversible contacts; 

- the combinatoric entropy due to the choice of which particular monomers (among those able to participate in 
bonds formation) actually do create bonds; 

- the loss of conformational entropy of the polymer chain due to bonds formation, and in particular, the entropic 
penalty of loop creation between two bonded monomers. 

Among a variety of macromolecular systems with thermoreversible interactions we pay a special attention to a 
class of so-called "RNA-like" polymers. These polymers arc distinguished from other biologically active associating 
polymers, such as, for instance, proteins, by a capability of forming only " clover leaf-like" (or "cactus-like") secondary 
structures. Indeed, the formation of a thermoreversible contact between two distant bonds in a RNA molecule (or 
in a single-stranded DNA) imposes a nonlocal constraint on a number of other possible bonds: all bonds in a RNA 
chain are known to be arranged in a way to allow only hierarchical cactus-like folded conformations topologically 
isomorphic to a tree. The pairs of bonds, which do not obey such a structure are called " pseudoknots" and are 
forbidden for RNA-type molecules. A quantitative definition of conformations allowed and forbidden in RNA-like 
polymers is described in the Section ^1 Note, that for proteins such a constraint is relaxed. 
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There is plenty of works concerning the statistical properties of associating polymers without any constraints on 
the topology of thermoreversible contacts 1^ . For such systems one can construct a satisfactory mean-field theory 
0. However RNA-like polymers are much less investigated. 

It is known that depending on the temperature (which comes into play via the temperature dependence of the 
statistical weights of the bonds) and on the primary structure of the chain, the RNA-like polymer may be found 
in the following four distinct states — see FigQ] First, at high enough temperature any RNA-like molecule forms a 
usual Gaussian or non-self-intersecting (depending on the quality of the solvent) coil with no, or almost no reversible 
bonds formed. Second, at lower temperature it may form a highly branched "genuine cactus-like" state topologically 
isomorphic to a randomly branched tree. Third, there is a strong evidence that at low enough temperatures the RNA 
molecules with highly irregular primary structure tend to form glassy states. Finally, in the case of some special 
primary structures, there exists an energy gap between the ground (so-called "native") energy state of the RNA 
and all the other states. The RNA can, therefore, in this case form a highly deterministic secondary structure (as 
compared to annealed random and quenched random states in the previous two cases). 
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FIG. 1: Schematic picture of possible thermodynamic states of RNA-like heteropolymers. 



In particular, if the primary structure of a RNA consists of two (or several) complimentary parts, it tends to form a 
double-folded (or, in case of several parts - a cross-like, as in typical tRNA) structure. Such a structure, contrary to 
a randomly branched one, is highly deterministic — one can predict easily which particular monomers will aggregate 
with each other — and has a different asymptotic (for large chain lengths N) mean-square end-to-end distance, R. 
For example, if one neglects the volume interactions (i.e. the non-saturating interactions between monomers which 
are far from each other along the chain), R^c ~ N^^^ for necklace (or "double folded") conformation, rather than 
^rb ~ iV^/'* in the randomly branching case. 

Let us give a sketch of what is known about the transitions between the four states depicted in the Fig^ First 
of all, the transition between the single-stranded and the double-stranded RNA (or DNA) commonly known as a 
helix-coil transition, which is a typical representative of the transition between the coil and the deterministic ground 
state, is rather well understood |7|, the modern developments in this field can be found in It is known that the 
particular form of loop penalties plays crucial role in whether it is a genuine phase transition or just a collective one. 
The contribution to the polymer partition function emerging from the loops is usually supposed to have the scaling 
form N"" (where N is the loop length) — see, for example, 7J. If < a < 1 (a = corresponds to the absence of any 
loop contribution) the passage from a coil to a helix is not a phase transition, but a cooperative phenomenon with 
a more or less (depending on the value of a and on the cooperativity of bond formation) sharp crossover. On the 
contrary, if a > 1 there is a genuine phase transition. The situation concerning randomly branched RNA molecules is 
more complicated. Despite the denaturation of a cactus-like RNA has been investigated for several decades, starting 
from classical works 0:13 1 '^^^ shown only recently ^3 that the character of transition in RNA resembles in some 
aspects the coil-helix transition in DNA. Namely, for large loop penalties (a > 2 for RNA) one has a genuine phase 
transition, while for smaller loop penalties there is just a crossover behavior. 

Note, that in both aforementioned cases the heteropolymer nature of real RNA structure played a minor role 
in the transition under discussion: one could have considered all monomers as identical and still have a proper 
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insight into what goes on. On the contrary, in the problem of glass transition in RNA-like polymers, the frozen 
heteropolymer structure of a chain plays a crucial role: the presence of frustrations in the system is itself a peculiarity 
of heterogenous systems There is, hence, a barest necessity for a statistical theory of heteropolymers forming 

RNA-like structures. The investigation of some thermodynamic properties of random heteropolymers is addressed 
in a few recent theoretical papers [TH ITsL H^. On the other hand, the theoretical results are usually obtained 
not for "real" heteropolymer with finite number of link types, but for polymers with frozen Gaussian distribution of 
association constants, sometimes even excluding the loop factor. It is believed (see, for example, jllj) that the loop 
penalties are less significant for understanding the physical properties of RNA-like molecules at low temperature (i.e., 
close to the glass transition) when the typical sizes of loops become rather small. 

We would like to stress that for the best of our knowledge, the theoretical consideration of the possible transition 
between regular (for example, double-stranded) and randomly branched cactus-like states of RNA is absent. It seems 
to be instructive, therefore, to consider a model of a RNA-like polymer which shares the peculiarity of a chain non- 
uniformity in a way, which can give rise to the formation of linear double-stranded structure, with an advantage of 
being exactly solvable. Presenting here such a model, we believe that besides of solving specific model of the transition 
between double-stranded and randomly-branched states of RNA, we also provide a better insight into the structure of 
possible states of heteropolymer RNAs, thus contributing in the construction of the general theory of heteropolymers 
forming RNA-like structures. 

m times n times 

Specifically, we consider in this paper a A„jBri = AA...A BB...B diblock copolymer, whose links can form local 
reversible bonds with each other. We show that, depending on the relative strengths of A-A, A-B and B-B bonds, such 
a polymer can be in one of two different states. Namely, if a self-association is preferable (i.e., A-A and B-B bonds 
are comparatively stronger than A-B ones) then the polymer forms a typical randomly branched cloverleaf structure. 
On the contrary, if alternating association is preferable (i.e. A-B bonds are stronger than A-A and B-B contacts) 
then the polymer tends to form a generally linear necklace structure somewhat reminiscent of a double-stranded DNA 
(with, probably, some rear side branches and loops, which, however do not influence the overall characteristics of the 
chain). In this paper we consider the transition between these two states (cloverleaf and necklace) in detail, and show 
that it turns out to be a 2nd order phase transition. 

The paper is organized as follows. In the Scction^we discuss the model under consideration, and derive the basic 
equations for the partition function. In Section fllll which plays the central role in the paper, we solve these equations 
in the approximation of no loop penalties, and study the case of large equal lengths of blocks, m = n ^ 1. We show 
that a necklace-cloverleaf transition is a 2nd order phase transition, compute the corresponding phase diagram and 
discuss the behavior of the free energy near the transition point. In the Section llVI we solve the basic equations with 
"ideal" loop weights {a = 3/2) and show that the presence of these loop factors does not have a strong influence on 
the phase transition. Finally, in the Section we summarize the results and briefly discuss the perspectives. 



II. THE MODEL 

Consider a A„Bm diblock copolymer whose monomers A and B are capable of forming thermoreversible bonds of all 
possible types A-A, B-B, and A-B. Let the statistical weights of bonds be equal ui, V2, and u, respectively. Having 
the applications to RNA molecules in mind, assume also that the structures formed by these thermoreversible bonds 
are always of a tree-like type, as shown in Fig|5^. It means that we restrict ourselves to the situation in which the 
chain conformations with " pseudoknots" shown in Fig|2l3 are prohibited. 

This definition can be formalized as follows. Take a chain C and enumerate its links: l,2,...,n. Consider two 
different thermoreversible bonds b(ii, 12) and &(ji, J2), where b{ii, 22) corresponds to a contact between two monomers 
ii and 12 (where 1 < «i < J2 ^ n), and &(ji, ^2) is formed by monomers ji and j2 (where 1 < ji < J2 < n). We 
say that the chain C forms a RNA-like structure if and only if for any two different bonds 6(ji,?2) and 6(ji, J2) the 
following inequality holds: 

(ji - ii){.h - - «2)(j2 - ^2) > (1) 

To make the difference between allowed and not allowed structures geometrically more transparent, redraw the 
structures shown in Fig|2jL and Fig[5j3 in the following way. Represent a polymer under consideration as a straight 
line with active monomers situated along it in the natural order, and depict the thermoreversible bonds by dashed 
arcs connecting the corresponding monomers. Now, the absence of pseudoknots simply means the that the arcs do 
not intersect — see Fig|21;,d. 
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FIG. 2: (a,b): Schematic picture of allowed (a) cactus-like and prohibited (b) pseudoknot configurations of the bonds; (c,d): 
Arc diagrams corresponding to configurations (a) and (b), respectively (note the intersection of arcs in (d); (e): The height 
diagram corresponding to the bond configuration (a). 

Moreover, we assume for simplicity, that except pseudoknots, all other bond configurations are allowed. This 
means, in particular, that we do not require any minimal loop length, as well as we do not yet take into account the 
cooperativity effect (the fact that if two links are connected with each other, then the two adjacent links have larger 
probability to be also connected). These assumptions are known to be false for real RNA molecules (for example, 
there are no loops shorter than 3 monomers in RNA chains ^3)- However, one can speculate that — see, for example, 
|15| — if the links of the chain are considered as renormalized quasi-monomers consisting of several "bare" monomer 
units, the made assumptions seem to be plausible. 

Let us also mention another possible representation of the secondary structures (configurations of bonds) allowed 
in the system under discussion. Consider again the arc diagram representation. Let then unconnected monomers, the 
monomers at which any arc starts and those where it ends correspond to vectors (1,0), (1,1), and (1,-1), respectively. 
Then any bond configuration is mapped to some sequence of vectors which can be considered as a random walk in 
(1+1)D (see Fig|2l3) situated in the upper half-plane. The walk begins and ends on the x-axis; the end points are 
separated by the distance N — the total length of the polymer. Such a representation is often called the "height 
diagram" of the secondary structure. Note, that the height of the point in the phase diagram is exactly equal to the 
number of arcs going above the corresponding point on the arc diagram, i.e. to the number of bonds one have to 
break to reach the corresponding monomer from the starting point of the chain. The critical exponent 7 connecting 
the mean height of such a diagram with the length of the polymer {{h) ~ N'^ , ^ 7 < 1), called roughness exponent, 
is an important characteristic of the state of the system. In the randomly branched homopolymer state 7br = 1/2, 
while in the glassy state of random RNA it was found numerically 7gi ~ 2/3, and in the recent work the 
analytical estimate 7gi ~ 5/8 has been obtained for the same system via the renormalization group procedure. 

Consider now the topology of the bonds formed in the AB-copolymer. One can easily see that A~B bonds separate 
the independent regions along the chain. Let us start moving along the chain from the left end and find the first A-B 
bond. As the arcs cannot intersect each other (see Fig|2), all the monomers under the given arc can be connected 
only with each other. Moreover, the A monomers in the beginning of the chain can form only A-A bonds (as we are 
speaking about the first A-B bond), and the B monomers in the rest of the chain can form only B-B bonds with each 
other, too. This separation of a chain into three independent parts allows us to write down the following Dyson-type 
equation for the partition function of the whole chain: 



G2{n,m,q\vi,V2,u) = G'(7i, q|vi) G(m, q|v2) +u 




(2) 



X G2 ° {n — i — l,m — j — l\vi,V2,u) 
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where G2(n, m, q|fi, V2, u) is the desired partition function (written in the q-space), G(fc, q|u) is the partition function 
of a self-associating chain of A or B monomers of length k (we assume for simplicity that the only difference between A 

and B monomers is in the weight of bond formation, respectively, vi and U2), and G2°°^\n, m\vi, V2, u) is a partition 
function of a diblock-copolymer which is forced to form a loop (i.e. its end monomers are forced to form a bond). 
Finally, the first term in the sum on the r.h.s. of eq. [3 is due to the fact that there can be no A-B bond at all. Note 
that to solve (0), apart from defining the homopolymer partition functions G one should allow for some particular 
connection between the partition functions of a "simple", G2, and "looped", G2°°^\ chains. 

Introducing the generating functions for all relevant partition functions: 

00 00 

G2{si,S2,q\vi,V2,u) = ^ ^G'2(m,n, q|wi,W2,u)s"s2; 

m— n— 

00 

G(s, q|«) = V G(m, q|«)s"; (3) 



00 00 



Gf°P^(si,S2ki,«2,u) = V VGf°P^(m,n|«i,t;2,u)« 



m— n—0 



we can rewrite equation |2Jl in the following form: 

G2(si, S2, q\vi,V2, u) = G(si, q|ui)G(s2, q|w2) + usiS2G{si,q\vi)G{s2,q\v2) 

xG^,'°°''''(si,S2|wi,t;2,ti) 



(4) 



(see Fig|3^ for the diagrammatic form of this equation). Moreover, the partition functions of a homopolymer forming 
clover-leaf structures (i.e., structures without pseudoknots) are known (see, for example, 0,0,^3) to satisfy the 
equation 

G(s, q|«) - g{.s, q) + vg{s, q)G(s, Ci\v)G^'°°^^ is\v) (5) 

(see FiginiD for its diagrammatic form) and the function g{s,q) is a partition function of a free chain (i.e., chain 
without any thermoreversible bonds or volume interactions). 
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FIG. 3: Diagrammatic form of Dyson-type equations: a) for Eq.Q, b) for Eq.(|SJ. Thin lines correspond to free propagators 
17(5, q), thick hatched lines — to the homopolymer propagators G{s,q\v), and thick grey lines — to the diblock copolymer 
propagators 6*2(81, S2, q|ui, i'2, w). Finally, open and filled circles correspond to the intra- and inter-species bonds, respectively. 

The particular form oi g{s,q) is known to be irrelevant for long enough chains (see, for example, 0). In our case, 
when a chain can form short-ranged loops, the particular form of the bare propagator infiuences importantly the 
entropy penalties for closing a chain into a loop. Note, however, that we introduce the loop penalties explicitly via 
the connection between G and G''°°P'' in (0 and therefore we assume the form of the function g{s, q) to be irrelevant. 
We can thus set without the loss of generality that 

where a is a mean length of a covalent bond between two subsequent monomers. After the Fourier-Laplace transform 
one gets 
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In what follows we chose the length units to set = 1. Note, that while the Gaussian form of the free chain 

propagator is a matter of assumption, the value of the propagator a.t q — (namely, (1 — s)~^) is independent of 
its particular spatial dependence. Indeed, the partition function of a free chain with given number of monomers 
Z{N) — J Z{N,r)(Pr is by definition equal to unity (see Q)- Therefore, the value of g{q,s) at g = 0, which is just 
Sw=o^(^)*^' equals always g{0,s) = 1/(1 — s). 

The crucial difference between equations (@J and ^ is that unlikely to the latter one, the former equation is linear. 
This corresponds to the aforementioned fact that the A-B bonds separate the chain into independent parts, therefore 
the global structure of the bonded chain is just a linear sequence of independent "blobs" separated by A-B bonds, 
only A-A and B~B bonds allowed inside each blob. This notion leads us to understanding of the main feature of 
the system under consideration: the possibility of a necklace-cactus transition. Indeed, for 2 2> u one can neglect 
the alternating A-B bonds and the chain forms just two independent "cactuses" consisting of A and B monomers, 
the roughness exponent 7 being equal to 1/2 similarly to the homopolymer case. On the other hand, if u ^ ui 2, 
the alternating bonds dominate in the structure, and the chain forms a "hairpin-like" structure with 7 = 1. Even a 
stronger conjecture is true: if the fraction of A-B bonds is nonzero in the limit of a long chain, the roughness exponent 
7 of such a configuration is equal to 1 . 

We are going now to analyze the described transition in more detail. To close the system of equations we should 
now address the aforementioned problem of the connection between "simple" and "looped" partition functions. In 
the two following Sections we consider, respectively, the following assumptions about this connection: there are no 
loop weights (Section lllip : there are "ideal" loop weights (Section llV|) . 

To be specific, in the Section IlIII we assume the most simple possible situation. Namely, we postulate: 

G^2°°^\si,S2) = G2(.Sl,S2,q 

G'(l°°P)(s) = 0(5, q = 0) 

Here the functions G(s,r) and G(si,S2,r) are the partition functions of the homopolymer and copolymer chains in 
the r-space. The assumption (jSJ corresponds obviously to omitting any entropic penalties for loop formation. We 
suppose that such an approximation can be valid at least for low enough temperatures, for the justification of this 
assertion see the beginning of the Section ITTll 

In the Section im we assume more natural connection between the partition functions under discussion: 

G'^'°°P^(si,S2) = j^-r^ I G2(si,S2,q)dq = G2(si,S2,r = 0) 
1 f 

G(i°°P)(s) = ^ j G(s,q)dq = G(s,r = 0) 

where the r.h.s. is written also in terms of the partition function in the r-space. Thus, we simply assume here that the 
partition function of a loop (i.e. chain, actually bonded at the end points) equals the partition function of a similar 
chain without end-to-end bond, but whose end monomers are situated at the same point in space. This assumption 
on the first glance looks almost obvious and it is hard to understand why it could fail. Note, however, that the 
above equations do not allow for volume interactions, i.e. interactions of non-saturating nature between monomers, 
situated far from each other along the chain, but close to each other in the real space. These interactions, however, do 
not influence the topology of possible cloverleaf conformations. Having in mind the discussion of the major factors, 
governing the the problem of a RNA-type folding, which we have presented in Section one can hope to take the 
volume interactions into account via altering the loop weight function, or, in our formalism, by choosing the proper 
connection between the "looped" and "non-looped" partition functions in the form different from the ideal one ©. 



III. NO LOOP WEIGHTS 

In this Section we calculate the partition function of the chain with no penalties for the loop formation. To justify 
this simplification let us remind that the loop weight is supposed to be relevant only if the loops formed in the system 
are long enough. Indeed, the loop weight, i.e. the probability that a polymer chain of length I forms a loop, is, at 
least in the first approximation, proportional to , were Ri is the typical end-to-end distance of such a chain, and 



0) = y G2(si,S2,r)dr 



(8) 



= / G(s, r)dr 
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D is the space dimensionality. In turn, Ri is proportional to some positive power of / depending on the state of the 
system. Therefore, the loop weights are of order of unity for short loops and much less than unity for long ones. In 
other words, the loop weights effectively suppress the long loops. If such loops are already suppressed (which is indeed 
the case if m, wi,2 3> 1, when almost all the monomers are involved into bonds), the particular form of loop weights 
is not much relevant. Hence, for large enough values of association constants (i.e., for low enough temperature) the 
no-loop-weights approximation seems to be realistic. 

Keeping this in mind, we can simplify equations l@J and jSJ as follows. We substitute the total partition functions 
G(q — 0) and G2(q — 0) instead of loop ones 

(^(loop) G^^°°P^ and then set q equal to zero. We arrive therefore to 

a set of two algebraic equations: 

G{s\v) = giO,s) + vg{0,syGHs\v) ^^^^ 

G2{S1,S2\V1,V2,U) = G{si\vi)G{s2\v2) + "Si S2G(si |wi )G'(s2 | ^2 )G2 (si , S2 | Wl , U2 , u) 

The solutions to read: 



G{s,v)^ i^— ^ 11 

for the homopolymer partition function (we have chosen the branch of solutions which approaches G = (1 — s)~^ as 
V 0); and 

1 

G2(si,S2) 



G ^isi,Vi)G-'^{S2,V2) - US1S2 

A ^ ^ > (12) 

AV1V2SIS2 




(1 - Si - V(l-Si)2-4t;is2) (1 - S2 - ^{l-S2Y-^V2i 

for the block-copolymer partition function. 

Note that though we have introduced the generating functions G(s) and G2(si, S2) as purely formal objects, one can 
ascribe a physical meaning to the result obtained in \i2\ . Indeed, imagine a polydisperse system of AB-copolymers 
in the equilibrium with a reservoir of A- and B-monomcrs. The length of A- and B-blocks will then be controlled 
by corresponding fugacities si and S2, and the partition function of such a block-copolymer is given by (O. An 
important note is that this situation would be fully annealed with respect to the lengths of the blocks, while in most 
of (at least biologically relevant) experimental cases the length and primary structure of a polymer are quenched. In 
what follows we are going to consider the quenched case when the length of the blocks is fixed and does not depend 
on the particular form of the interaction between monomers. 

It seems now instructive to study the analytic structure and determine singularities of H12|) . because the particular 
form of these singular points is responsible for the equilibrium conformations of the system under consideration. 
Indeed, one can easily see that there are three possible singularities of 112|l : two square-root singularities at the 
points 

-(^)^ -(^)^ 

and the simple pole singularity at the point 

4tiiW2SlS2 



_ _ ^(l-si)2-4«is2j [1-S2- ^{l - S2Y - ^V2si^ 



(14) 



The square-root singularity corresponds (see, for example, [Tfll) to the formation of a usual randomly branched 
clover-leaf structure with 7 = 1/2 due to the self-association of monomers of type A (or B), while the simple pole 
signals the formation of a linear (7 = 1) double-folded conformation due to the alternating A-B association. Thus, if 
the square-root singularity "wins" (i.e. becomes more important), the global conformation of a chain is a randomly 
branched, while if the pole singularity "wins", the global conformation is linear. One can imagine also a situation 
when (say, for a very long A-block and short B-block) there can arise also a mixed situation, when both singularities 
are of equal importance or, saying it in more physical words, when within a single polymer there is a phase equilibrium 
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between a necklace (double-stranded) phase and a homopolymer (in this case - A) cloverleaf phase. However, as we 
show below, this situation actually never exists in the particular system under discussion. 

If a system is controlled by one single activity s, the meaning of "winning singularity" is well-defined: in the 
thermodynamic limit the singularity closest to s = is dominant. In the case under consideration we have two 
independent activities si and S2, which makes the situation somewhat more tricky. However, one can handle this 
problem in the case of two activities considering the thermodynamic limit: m = aN; n = (1 — a)N when N oo 
and a = const. 

Here we present only the results, while the details of the algorithm we are using are explained in the Appendix A, 
and the details of particular computations - in the Appendix B. 

Depending on the particular value of the association constants {u,vi,V2} either the square-root or the pole singu- 
larity become dominant. So, we have two options: 

If u < Utr = y/viV2, then the square root singularity is dominant and the system is in the cloverleaf phase with the 
free energy per link of the chain / asymptotically (when — s- oo) equal to 



clov 



Ins! 



-aln(l + 2V^) - (1 - a) ln(l + 2y^) 



(15) 



(See Ea.(lB3ll^ 



If u > utr, then the pole singularity is dominant and the polymer is in the necklace conformation. The corresponding 
free energy reads 



/neck = lns;,i, = -aln(l + 2V^cosh(i? - 0) -(!-«) ln(l + 2V^cosh(i5 + 0) 



(see HB9|l'). where 



4 ViV2 



(16) 



(17) 



and ^ in (|16|l is the solution of the equation 

4V?r^sinh 2^ + 4(1 - 2a)Vwr^sinh 2i? - a/^sinh(i? - + (1 - a)^/v^smh{-d + ^) = (18) 
on the interval (— i?, i?). 

In Figure ^ we have plotted the free energy of the system defined by equations H15|l and (|16|l as a function of 
variables u, vi and a — see Figs. ^,b,c, respectively. 




FIG. 4: Free energy of the system defined by (1151 and HIGH as a function of u (a), vi (b) and a (c). (a, b): Thin lines correspond 
to vi = V2 = 1, thick lines - to i;2 = 2, vi — 0.5, dashed lines - to V2 = 5, vi = 0.2. (a) Curves corresponding to a = 0.1, 0.5, 
and 0.9 are marked by letters a, b, and c, respectively, (b) Curves corresponding to u = 1 (the transition point), 10, and 100 
are marked by letters a, b, and c, respectively, (c) Thick lines correspond to u = 2, V2 = 1, thin lines, both in the main figure 
and the inset - to u = 2, V2 = 20. Curves marked by 1, 2, 3, 4, and 5 correspond to a = 0.9, 0.7, 0.5, 0.3, and 0.1, respectively. 
Note that the thin lines intersect at the same point v\ = V2 which belongs to the cloverleaf phase. This is not true for the thick 
lines since the point vi = V2 belongs to the necklace phase in this case. 
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The free energy does not depend on u in the cloverleaf phase and it diverges quadratically at the transition point 
from the cloverleaf value, thus giving rise to the discontinuity of its second derivative. For /(wi) dependencies this 

discontinuity is less evident. To make it more clear we have plotted in the inset the function ^^-^^p^. The transition 
from the cloverleaf to the necklace phase is, thus, a 2nd order phase transition. Moreover, it is very important that 
the function /(a) is always a convex curve. This fact provides a support for our assertion that there is no any phase 
separation inside the system under consideration. Indeed, if there were the regions where the free energy as a function 
of a is concave, the system would have minimized its energy in a usual way, by separating into two necklace-phases 
with different effective values of a. However, as we see, it is not the case. Together with our consideration in the end 
of Appendix B where we prove that there is no necklace-cloverleaf phase coexistence in the system, this leads us to 
the conclusion that the homogenous phase in the system under consideration for any a is always more stable (i.e., 
has lower free energy) then two phases. 



IV. IDEAL LOOP WEIGHTS 



In this section we analyze the exact solution of Q and |(SJ| with the loop weights defined by ^ and with the 
propagator of a free chain given by Q) . Substituting Q) and ^ into Q one gets a nonlinear integral equation which 
is (see 0,^3) solvable by means of integration over q. One gets 

G{s,v,q)^ ^ (19) 

el —X[s,v) 

(compare to the "bare" function {Tj) and recall that we have set a^/6 = 1) where X{s, v) satisfies the following equation 



(4^)3/2x(A-s) =Vs2Li3/2(A) 



(20) 



and Li3/2(A) is a polylogarithm function of order 3/2. Equation H2()|l always has a solution for small enough s. When 
s is increasing, H2U|I loses the solution at the point defined by the derivative of H2U|) with respect to X: 



{At:)^/^X{2X - s) vs^UiMX) 



(21) 



Let us note that at this point the solution X{s,v) (and, therefore G{s,v,q)) has a square root singularity similar to 
that of equation H13|l . Substituting H19|) together with the first equation of ^ into one arrives at the partition 
function of a diblock copolymer with ideal loop weights: 



A2{si,S2) ^ j^-:^ y G2(si,S2,q)dq 



V" 



X{si,Vi) - X{S2,V2) 
X{si,Vl) - Si X{S2,V2)-S2 



US1S2 



visl 



V2SI 



(22) 



The equation (|22|l has a simple pole singularity at the point when the whole expression in brackets in the r.h.s. equals 
zero. Note that there is no singularity due to the denominator of the ratio inside the brackets, since the nominator 
and denominator always tend to zero simultaneously, while the ratio as a whole stays finite. Thus, one sees again 
that the conformation of the block-copolymer under discussion is governed by an interplay between a square root 
singularity (corresponding to the cloverleaf phase) and a simple pole (corresponding to the necklace phase) . 

Similarly to the previous section, we sketch here the results for the free energy and the transitions in the system, 
putting all the details of computations into Appendix C. 



The transition point now is given by 



U = 



V2S*2 {Xl{s\r^ - 1) - Visl (XI(S^)-I - 1) 



(23) 



where {s\,Xl) and (52;^!) f^e common solutions of equations (pn|l and for v = vi and W2, respectively. In 
Figl^lwe plot the phase transition line C/ti (i^i, ^'2) as a function of Vi/v2 for different values of the mean homopolymer 
association constant ^viV2. Note that in this case the dependence C/tr(i'i/i'2) is more rich than in the no- loop- weight 
case, where the transition point was at Utt = = 1 independently of the ratio of the homopolymer association 
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constants. In particular, this nontrivial dependence means that there is a possibihty of reentrant cloverleaf-necklace- 
cloverleaf transition in the system under discussion. Indeed, it is natural to assume (compare to 0,0, UJi) that the 
association constants depend on temperature in a usual Boltzmann way: 

u = uoexp{-eo/{kBT)); Wi,2 = «d' ^ exp(-ei^2/(fcBr)) (24) 

This means that the change of temperature in any particular system with given chemical structure and given values 
of uq, Uq'^, £0,1,2 corresponds to a movement along some straight line on the log-log plot in Fig|S| We easily see that 
such a straight line can easily have at least two intersections with the phase transition line thus giving rise to the 
aforementioned reentrant transition |2J| . 




FIG. 5: Phase diagram of the cloverleaf-necklace phase transition with ideal loop weights. Thin, dashed and thick lines 
correspond to the product viV2 = 100, 10, and 1, respectively. 



Moreover, the free energy of the system with ideal loop weights reads 

/ = aln4 + (l-a)ln4 (25) 

where in the clover leaf phase the acting singularity values a-i'e equal to the 2 defined in (|23|l . while in the 
necklace phase they are the solutions of the set of equations (|26|) , (jCip , and (jC7p lying within the rectangle < < 

s*i;0< 4 < si 

In Figl^we plot this free energy / as a function of a for several different values of wi, V2 and u. One can notice the 
similarity between this plot and Fig^jD. The curves are once again always convex, so there is once again no internal 
phase separation in the system. Moreover, one can easily check that the necklace-cloverleaf transition is again of 
2nd order. Thus we see that the difference between no-loop-weight and ideal-loop-weight cases is mainly qualitative, 
except the peculiarities in the behavior of the transition point Utr{vi, V2) and the possibility of the reentrant transition 
in the latter case. 



V. CONCLUSIONS 



We have shown that a diblock copolymer capable of forming RNA-like structures due to reversible bond formation 
can exhibit a phase transition between two different thermodynamically equilibrium states, a "necklace" with the 
roughness exponent 7nc = 1, and a "cloverleaf with 7i-b = 1/2. The transition is governed by the values of equilibrium 
association constants, and the weights of loop formation in the system. 

We have demonstrated that a cloverleaf-necklace transition in the thermodynamic limit is a 2nd order phase 
transition both in the case of no- and ideal loop weights (in these cases we have calculated the free energy of the 
system exactly). One can expect that allowing larger penalties for loops formation, the new phase, that of a coil, may 
arise in the system, similarly to what predicted in jlO'l, and the order of the necklace-cloverleaf transition itself may 
change. 

The interesting feature of the systems under discussion in the case of ideal loop weights is the possibility of a reentrant 
necklace-cloverleaf transition for some particular values of parameters determining the temperature dependence of 
the association constants. 
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FIG. 6; Free energy of the system with ideal loop weights as a function of a. Thin, thick, and dashed lines correspond to 
vi — V2 — 14.1; vi = 10, V2 — 20; and vi = 4, V2 — 50, respectively. The curves corresponding to u — 40, u — 80, and u = 400 
are marked by a, b, and c, respectively. Note that u = 40 (curves a) correspond to a point in the necklace phase just above the 
transition. In the cloverleaf phase (see IIU2t 'l the free energy is linear in a. 



Moreover, the proposed approach allows to calculate the mean free energy of a frozen ensemble of the diblock 
RNA-like copolymers with some given distribution of bond lengths. 

Finally, we would like to emphasize that the method of calculating the partition function of a diblock-copolymer 
based on iQJ , Q (or equations (jTUIl in the case of no loop weights) can be straightforwardly generalized to the case of 
alternating polyblock RNA-like copolymers of type Af^Bf^ Af3B4^...At^_^Bf^ with all lengths of blocks being different. 
This situation models a particular example of an alternating heteropolymer RNA-like molecule with quenched primary 
structure defined by the sequence of length {ti, t2, t^}. For example, the case of a 4-block copolymers seems to be 
of particular interest, as in this case one can expect the equilibrium cross-like conformations reminiscent of that of 
the tRNA. 

The corresponding general Dyson-type equation for the generation function G„(si, S2, Sn), where 

oc oo oc 
ti=0t2=0 tk=0 

takes the following form (compare to @): 

fe 

Gfc(si, Sfc) = Gi(si) Gfc-l(s2, Sfc) + ^ Si Sr Gi(si) Gr(si, S^) Gfc+l-r(Sr, Sfc) (26) 

r=2 

which has the diagrammatic structure shown in FigEI (compare to Fig|3Ji). The equation |SJ|, as well as it diagrammatic 
structure shown in Fig|2|D remain without changes. 



r 




FIG. 7: Diagrammatic structure of equation 1261 . 



We believe that equation (|26|l can be served as a basis of a renormalization group procedure, thus connecting the 
exact partition function, G(4), of an 4-block copolymer with an approximative expression for G(4) constructed of 
two diblock copolymers, G(2), i.e. G(4) « G(2) x G(2). It is obvious that if all couphng constants are identical, then 
we have a RNA-like homopolymer and the the last relation is exact. However if A-A, B-B and A-B constants are 
slightly different, we should have G(4) — G(2) x G(2) = £{u, Wi, ■^2) 7^ 0. Considering e as a perturbation we may look 
how e is transformed under the standard coarse-graining procedure. This work is in progress and the corresponding 
results will be reported elsewhere [2l| . 
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APPENDIX A: SINGULARITIES OF THE GENERATING FUNCTION G{si,S2) 



In this Appendix we address the procedure of calculating the asymptotic behavior of the coefficients of a generating 
hmction (compare to 

oo 

G(si,S2)= J2 (Al) 

n,m— 

in the limit m + n = N ^ oo and m/{m + n) = a = const. The procedure suggested below is similar to that 
developed in and we advise the reader to address this reference for the accurate mathematical proofs. However, 
the algorithm under discussion is not of a common use, so it seems instructive to give a short explanation on the 
physical level of rigor. Note, that we pay attention only to the exponentially increasing multipliers in the asymptotic 
behavior of the free energy, thus neglecting the logarithmic corrections. 



First of all, let G(si, S2) defined by (|Aip be an analytic function of si, S2 in the circle ■\/|sip + |s2p < e for some 
e > near the origin si = S2 = 0. Assume also, as it is usual in the statistical mechanics, the dominant singularities 
of G with respect to each of the variables (for fixed another one) to be real and positive. 

Now, it is known that to calculate a function Gl{s) defined as follows: 

oc 

GJ(s) = ^aMs2' (A2) 

whose dominant singularity defines the asymptotic of ai^i at i 00, one can use the following trick (see, for example, 
[2^). Consider G (sz, |) as a function of a variable z. If s is small enough this function is analytic with respect to z 
inside some ring where it can be represented by the Laurent series 



G 



00 

(sz,^) = J2 (A3) 



fc— — 00 



Note now that the coefficient feo(s) in this expansion is exactly equal to the desired function Gl{s) and therefore this 
function can be evaluated using the Cauchy theorem: 

Glis)^l-fG(.szrl) ^ (A4) 
Zm J (J \ z) z 

where the curve C encloses the point z = exactly once and belongs to the region of analyticity of the function 
G(sz,f). 

This consideration can be easily generalized in the following way. Suppose, we are interested in calculating the 
function 



where p and q are two mutually simple natural numbers. To perform the summation in (jA5|) . we consider a function 
G(sz*, 32"^*) and take the coefficient of order in its Laurent series: 
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where the integration curve belongs again to the region of analyticity of G{sz'^, sz p). 

This formaUsm allows one in principle to calculate the generating function G*(s, ^) for any given ratio of block 
lengths of the copolymer under discussion. However, on the one hand, the explicit calculation of G* is still rather 
cumbersome, and on the other hand, it is actually not needed. What we are really interested in is the asymptotic 
behavior of the coefficients of G*, or, in other words, the position and the nature of its dominant singularity (or 
singularities). The task of finding these singularities can be accomplished without explicit knowledge of G*. 

To do that, let us first introduce some additional definitions. Let H be the 3D surface of singular points of G (it may 
consist of several branches, of course) in the 4-dimensional space (Resi,Res2,Imsi,Ims2) and h be the intersection 
of H with the "real positive plane" (i.e., the quarter-plane of (si, S2) where Arg(si) = and Arg(s2) = 0). Now, the 
surface H separates the space (si, S2) into several domains. We pay a particular attention to the shape of the domain 
containing the origin (si = S2 — 0). Let us call its edge Hi {Hi C H). It follows immediately from the condition of 
single-variable dominant singularities to be real and positive, that for any given absolute value |si|, |s2|, the point of 
the edge of the domain Hi closest to the origin is the point with Arg(si) = Arg(s2) = 0, and, thus, it belongs to the 
curve h. In what follows we can, restrict ourselves to consideration of real and positive singularities of G. 

In a (si,S2) quarter-plane there exists some domain of analyticity of G which includes the origin, its edge is 
designated by hi. An example of such a domain of analyticity of G in a quarter-plane (si, S2) including the origin, 
and the edge, hi {hi C Hi), of this domain are shown in Fig|Hl 




FIG. 8; Examples of possible domains of analyticity. The curves ljA7|l corresponding to different values of a {ai < a2) are also 
shown. The touching points may (b) or may not (a) depend on a. 

Let us fix now some "direction" a — We are interested in calculating the generating function G*{s,a) 

corresponding to this particular direction as defined by ljA5|l . We are thus to integrate the function G{sz'^ , sz^p) 
as prescribed in HA6|I . By fixing some particular value of s, we define therefore a connection between the values of 
si = sz'' and S2 = sz^^ in the integrand as follows: 

= (sPzP'')(s«z-P«) = (ors?4-"=.s) (A7) 

Obviously, the corresponding line in the (si, S2) plane always intersects our domain of analyticity if s is small enough. 
Therefore, for small s there is no problem to choose the appropriate circle of integration \z\ — const. One should 
just choose 1^1 in such a way that the point (sj^j'', sjzl"'') belongs to the domain of analyticity. However, when one 
increases s, the intersection between this domain and the line ljA7|l may disappear at some point. In terms of the z 
plain, where the integration prescribed in (|A6I) is to be accomplished, it means that two singularities of the integrand, 
one being inside the integration curve, and another - outside, approach each other, and finally merge. Note, that 
they do merge indeed rather then just have equal absolute values since we have assumed these singularities to be real 
and positive. Thus, we see that the integral in ljA6|l has a singularity exactly at this point. 

Thus, summarizing the said above, we suggest the following procedure to find the singularities of the function 
G*{s,a): 

1. Find the lines of singular points h of the initial function G(si, S2) in the real plane; 

2. Define the domain of analyticity D which includes the origin si — S2 ^ and find its edge hi; 

3. For each value of < a < 1 find the value of s for which the intersection of D and the curve (|A7p ceases to 
exist for the first time. If hi is a smooth curve it means just finding a common tangent line of hi and ljA7|l . 
if there is a finite number of peculiar points on it (say, kinks or intersections of the branches), it includes also 
finding the conditions that (|A7I) contains these particular points. 
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APPENDIX B: CALCULATION OF THE FREE ENERGY IN THE NO-LOOP- WEIGHTS LIMIT. 



In this Section we apply the strategy presented in the Appendix A for calculating the free energy corresponding to 
the grand partition function (|12|l . i.e. that of the diblock copolymer with no loop weights. 

As prescribed by the suggested procedure, we should first find the domain of analyticity which includes the point 
si = S2 = 0. In Figiniwe have drawn these regions in the (si, S2) plane for some particular sets of vi, V2 and different 
values of u. The edges of the region correspond to singularity lines defined by H13|l and (|14|l . The square root 
singularities correspond to some straight lines si = \s^^'^\ and S2 = \s^'^\ (where one should choose the solution of H13() 
with the lowest absolute value), so for any (si, S2) with absolute values the rectangle (0 < |si| < js^'^'^l, < |s2| < I) 
the square roots in (|12|l are both analytical functions of their respective variables. Moreover, the pole singularity 1)14(1 
defines some other curve in the (si,S2) plane. This curve can cither have no branch within the rectangle described 
above, and than this rectangle is the desired domain of analyticity itself, or it can intersect the rectangle (as shown 
by lines a-e in FigO and then one should take the part below and to the left of it as a desired domain of analyticity. 



0.2 T 




FIG. 9: Examples of the analyticity domains of 1)12^ for «i — 1,V2 = 5. Rectangular domain with borders Si,2 = (1 + 2y«i,2) ^ 
corresponds to values of u below the necklace cloverleaf-transition (i.e., u < ittr = V5). The curves marked by a, b, c, d, e are 
the pole lines corresponding to u = 5, 10, 20, 50, and 200, respectively. 



The further calculations can be much simplified by introducing the new variables 

■51,2 



ti,2 = 



1 - Sl,2 



The square root singularities now correspond just to ti,2 = 1. The equation for the pole singularity 1(14(1 can be 
rewritten as follows: 

U=^^ = (Bl) 

(i-yr^)(i- v/r^) 

Now, since the function 

m ' 



decreases monotonously from +00 to 1 with t increasing from to 1, it is obvious that the solution of ((Bip in the 
square 0<ii<l, 0<i2<l exists for and only for [/ > 1, i.e. for u > utr — -Jv\V2. Moreover, we easily see 
that for u < uti when the region of analyticity is just the rectangle shown in FiglHl the corridor of analyticity for the 
integrand in the integral expression for G* l|A6p disappears exactly at the point when the line l(A7p passes through 
the point (s^^^, s^^^): 

The free energy per link of the chain is given therefore by 

/ciov = lns:q, = -aln(l + 2^1) - (1 - a) ln(l + 2^) (B3) 
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and the singularity of G* is in this case of square root nature which justifies our assertion that for u < utr the polymer 
is in the cloverleaf state. 

In the case when u > un the region of analyticity of ljA6|l disappears when the curves (|A7p and (|B1|I have a common 
tangent line. To proceed further let us introduce the new variables 4>i^2 as follows 

sin (/>i 

(pi = arcsin^i; Si = — , , i ^ 1,2. (B4) 

sm0i + 2y/vl 

The equations (|A7|I and (|B1(I can be now rewritten as follows 

C/"^ tany tany; (B5) 

and 

s = (l + 2V^/sin,/)i)-"(l + 2V^^/sin02)""' (B6) 

were (f)i^2 € [0, n/2]. 

The condition to have the common tangent line gives rise to the additional equation 

sin0i 1 — a /^sin(/>i cos02(sin0i + 2y^) 



sin 02 a \j V2 sin (j)2 cos 4>i (sin 02 + 2y/v^) 



(B7) 



One can easily exclude two of the three unknown variables (s, (j>i, (j>2) from equations ljB6p - ljB7|) . and arrive finally 
to the following equation 

4VwI^sinh(20 + 4(1 - 2a) Vwr^sinh(2i9) - a^smh{d - + (1 - a)V^sinh(i9 + (B8) 

for a new variable ^ = — ^InC/ — In tan with = — |lnC/ and additional condition — z9 It is easy to 

check that for any d > (i.e., for any u > un) the left hand side of this equation is an increasing function of ^ 
and changes the sign inside the interval {—'&,{}). The equation (|B8I) has therefore exactly one solution with respect 
to ^. The fact that this solution always stays within the interval (— i?,!?) and never riches its limits means that the 
dominant singularity in the necklace phase is always (independently of u and a) a regular point of the pole line, not 
the intersection of the pole and square-root line. This means that there is no coexistence of the necklace and cloverleaf 
phases in the system for any u and a. Moreover, after substituting this solution ^* into the original equation for the 
singularity of s (jB6|l we obtain the final expression for the free energy in this (necklace) phase: 

/nock = Indole = -aln(l + 2V^cosh(i? - C)) -(!-") ln(l + 2V^cosh(i? + C)) (B9) 



APPENDIX C: CALCULATION OF THE FREE ENERGY IN THE IDEAL-LOOP-WEIGHTS LIMIT 



In this section we apply the procedure suggested in the Appendix A to the calculation of the free energy of the 
block-copolymer with "ideal" loop weights. To that end we investigate the singularities of the grand partition function 
with ATi, ^2 defined by equation if^ . 

As we have mentioned already, there are two possible types of singularities of this partition function. The first type 
corresponds to the singularities of the homopolymer partition function, is of square-root nature and corresponds to 
solution of equations H2l)|l and (|21|l . The second one is a pole singularity corresponding to the zero of the denominator 
in the integrand on the r.h.s. of (|22|l : 

(Xi - X2) f - ^^"l '-us,S2 = (CI) 

V Visf V2S2 J 

where Si and Xi are connected via the equation 12U|) . Contrary to the first one, this singularity is of collective nature, 
i.e. it depends simultaneously on the values of both si and S2, and it also depends on the value of u. 

Moreover, in full analogy with the no-loop-weights case, the square-root singularities define for small values of u 
the region of analyticity of a rectangular shape in the (si,S2) plane. In this case (i.e. below the necklace-cloverleaf 
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transition) the acting singularity point is, independently of a, just the intersection point of the two square-root 
singularity lines, and the free energy equals to 

/ciov(wi,2, a) = -alnsl - (1 - a) Ins^ (C2) 

(compare to (|B3j) where si 2 are the solutions of H20() and H21|l with v = vi,V2- 

Now, with increasing u the line of poles enters the rectangle of analyticity. It first appears at the intersection point 
of two square-root singularity lines, thus giving rise to the following equation for the necklace-cloverleaf transition 
point: 



(^r - ^2) 



(C3) 



where we have inserted the solutions of the set (|2UI) . 121|l into (jClll . This equation leads after some algebra to (|23|l . 

To investigate properly the phase above the transition point (i.e., for u > Wtr) we have: i) to find the common 
tangent lines of the pole line (jCl|l and the line IjATp . and ii) to check whether there is a possibility for the acting 
singularity to lie at the intersection of the pole line and the square-root singularity line, i.e. to check whether there is 
a kink at this intersection. To do that we should first of all find the derivative of the pole line. The direct calculation 
gives: 



dsi 



x'Ai 



US2 
VlSl 



US2 



X. 



S2 



V2S\ 



visl 



= ds2 



x'2 



1 



USl 

V2S2 



USl 



X, 



v\s\ 



X2 

V2SI 



where X( 2 are the derivatives of Xi^2 with respect to si^2'- 

dX,, „ (47r)3/2x, + 2w,s,Li3/2(X,) 



X' = 



= x,. 



i = 1,2. 



(C4) 



(C5) 



ds, (4^)3/2X,(2X, - Si) - V^s|Uy2{Xi, 

and the auxiliary conditions (|20|l are taken into account. 

Now, let us notice that along the square-root singularity line (say, for si = s^), the corresponding derivative {X[ in 
this case) diverges, while, if we are not approaching the point (s*, S2), all the other coefficients in ljC4(l remain finite. 
We conclude immediately, that to comply this equation, the derivative ^ along the pole line should also diverge in 
the vicinity of the point si = Si. Thus, the pole line has a vertical asymptote at si = s^ and there is no kink at 
the intersection of the pole and square root lines (the case S2 — S2 absolutely similar). We see therefore that the 
acting singularity never lies at this intersection. This proves, similarly to the no-loop-weight case, that there is no 
possibility of the intramolecular necklace-cloverleaf phase separation in the system. 

To calculate the free energy of the necklace phase as a function of the composition variable a we should now equalize 
the pole line derivative (|C4|) and the derivative of the line ljA7p which is given by 



dsi 
ds2 



a si 

' S2 



(C6) 



Collecting equations 1C4|I and (|C6(I together, we get 



as2 



X'2 



USl 

V2S2 



USl 



Xi-si 
visf 



X2 

V2S\ 



(a - \)si 



X'i[\ 



US2 
VlSi 



US2 



X2 



S2 



V2S\ 



Xl 
vis\ 



(C7) 



This equation, together with H20() and lCip completes the set which defines four unknown variables si, S2, ^1, X2 at 
the pole singularity point. 

The free energy is again 

/neck = -a In si - (1 - a) In S2 (C8) 
It is calculated numerically and the corresponding plots are shown in FigElat length of the Section Hvl 
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